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20.  ABSTRACT 

General  Plasma  Theory  and  Simulation 

A.  '  The  Pierce  diode  linear  behavior  with  external  R,  C,  or  L  was  verified  very  accurately  by  particle 

simulation.  Final  report  issued. 

B.  The  Pierce  diode  non-linear  equilibria  with  R,  C,  or  L  are  described  theoretically  and  explored 
via  computer  simulation.  Final  report  issued. 

C.  ~  A  simple  model  of  the  sheath  outside  the  separatrix  of  an  FRC  was  modeled  electrostatically  in 

2d  and  large  potentials  due  to  the  magnetic  well  and  peak  were  found;  these  may  explain  the 
anomolously  high  ion  confinement  in  the  FRC  edge  layer.  Memorandum  to  be  issued. 

D.  -  A  planar  plasma  source  with  cold  ions  and  warm  electrons  produces  a  source  sheath  with 

sufficient  potential  drop  to  accelerate  ions  to  sound  velocity,  which  obviates  the  need  for  a  Bohm 
pre-collector-sheath  electric  field. 

Plasma-Wall  Physics,  Theory  and  Simulation 

A.  '  Collector  Sheath,  Presheath,  and  Source  Sheath  in  a  Collisionless,  Finite  Ion  Temperature  Plasma; 

final  report  in  preparation. 

B.  Potential  Drop  and  Transport  in  a  Bounded  Plasma  with  Ion  Relfection  at  the  Collector;  final 
report  in  preparation. 

C.  Potential  Drop  and  Transport  in  a  Bounded  Plasma  with  Secondary  Electron  Emission  at  the  Col¬ 
lector;  final  report  in  preparation. 

D.  A  movie  has  been  made  displaying  the  long-lived  vortices  resulting  from  the  Kelvin  ,  .elmhollz 
instability  in  a  magnetized  sheath. 

Code  Development  and  Software  Distribution 

A.  A  relativistic  Monte  Carlo  binary  (Coulomb)  collision  model  has  been  developed  and  tested  for 
inclusion  into  the  electrostatic  particle  simulation  code  TESS. 

B.  Two  direct  implicit  time  integration  schemes  are  tested  for  self-healing  and  self-cooling  and 
regions  of  neither  are  found  as  a  function  of  At  and  A*  for  the  model  of  a  freely  expanding 
plasma  slab. 
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SECTION  I:  GENERAL  PLASMA  THEORY  AND  SIMULATION 


A.  The  Pierce  Diode  with  an  External  Circuit.  I.  Simulations  in  the  Linear  Regime 
William  S  Lawson 

Abstract 

The  linear  theory  for  the  eigen  frequencies  of  the  extended  Pierce  diode  for  resistive,  capacitive, 
and  inductive  external  circuits  is  tested  via  computer  simulation  with  the  PDW1  code,  and  is  verified  to 
within  the  accuracy  of  the  simulation  ('between  I %  and  .01 7c  depending  on  the  circumstances).  The 
simulations  were  remarkably  inexpensive  in  terms  of  computer  time,  and  could  have  been  performed  on 
a  personal  computer. 

This  work  is  to  be  issued  as  an  ERL  Memorandum. 

B.  The  Pierce  Diode  with  an  External  Circuit.  II.  Non-Linear  Equilibria 

William  S.  Lawson 

This  work  is  to  be  issued  as  an  ERL  Memorandum.  The  abstract  follows. 

C.  Electrostatic  Potential  Formation  Due  to  a  Large  Dip  in  the 
Magnetic  Field  with  Application  to  FRC  Confinement 

S.  E.  Parker 

This  work  is  to  be  issued  as  an  ERL  Memorandum.  The  introduction  follows. 

D  Ion  Acceleration  in  the  Plasma  Source  Sheath 
C.  K  Birds  all 


B  .  THE  PIERCE  DIODE  WITH  AN  EXTERNAL  CIRCUIT.  II. 


NON-LINEAR  EQUILIBRIA 


William  S.  Lawson 
Plasma  Theory  and  Simulation  Group 
Electronics  Research  Laboratory 
University  of  California 
Berkeley,  CA  94720 


ABSTRACT 


The  non-umfcrm  equilibria  of  the  classical  (short  circuit)  Pierce  diode  and  the 
extended  (series  RLC  external  circuit)  Pierce  diode  are  described  theoretically  and  explored 
via  computer  simulation.  It  is  found  that  most  equilibria  are  correctly  predicted  by 
theory,  but  that  the  continuous  set  of  equilibria  of  the  classical  Pierce  diode  at  a=2rr  are 
not  observed.  The  stability  characteristics  of  the  non-uniform  equilibria  are  also  worked 
out.  and  are  consistent  with  the  simultations. 
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Electrostatic  Potential  Formation 


Due  to  a  Large  Dip  in  the  Magnetic  Field 


with  Application  to  FRC  Confinement 


S.  E.  Parker 

Plasma  Theory  and  Simulation  Group 
Electronics  Research  Laboratory 
University  of  California 
Berkeley,  CA  94720 


ABSTRACT 


Electrostatic  effects  may  explain  the  enhanced  plasma  confinement  observed  in  the  open  field  line 
region  of  field-reversed  configurations.  We  are  studying  the  formation  of  electrostatic  potentials  due  to 
nonuniform  magnetic  fields.  This  type  of  phenomenon  is  being  modeled  using  particle  simulation.  Our 
two  dimensional  (x,y)  electrostatic  particle  code,  ES2  has  been  used  to  study  the  formation  of  electros¬ 
tatic  potentials  due  to  magnetic  wells  and  peaks  with  an  injected  flux  of  well  magnetized  electrons  and 
weakly  magnetized  ions.  These  results  are  being  used  to  help  interpret  the  electrostatic  potential  con¬ 
tours  outside  the  separatrix  of  a  FRC,  which  are  really  cylindrical,  (r,z).  These  potentials  may  explain 
the  anomalously  high  ion  confinement  in  the  FRC  edge  layer.  Our  work  naturally  follows  from  the 
qualitative  description  of  electrostatic  effects  on  FRC  confinement  given  by  L.  C.  Steinhauer. 


D .  Ion  Acceleration  in  the  Plasma  Source  Sheath 
C.K.  Birdsall 

This  note  is  a  calculation  of  the  potential  drop  for  a  planar  plasma  source,  across  the  source 
sheath,  into  a  uniform  plasma  region  defined  by  E_=0  and  or  perhaps  d2<t>/dx2  =  0.  The  calculation 
complements  that  of  Bohm  (1949)  who  obtained  the  potential  drop  at  the  other  end  of  a  plasma,  at  a 
planar  collector  sheath.  The  result  is  a  relation  between  the  source  ion  flux  and  the  source  sheath 
potential  drop  and  the  accompanying  ion  acceleration.  This  planar  source  sheath  ion  acceleration 
mechanism  (or  that  from  a  distributed  source)  can  provide  the  pre-collector-sheath  ion  acceleration  as 
found  necessary  by  Bohm. 

[The  analysis  is  elementary  and  may  be  in  early  works,  such  as  Tonks  and  Langmutr  (1929)  and 
Harrison  and  Thompson  (1959).  Our  apologies  to  such  authors.] 

MODEL.  POISSON  EQUATION 

The  model  has  a  planar  source  of  both  warm  electrons  and  cold  ions  shown  in  Fig.  1.  Our  object 
is  to  obtain  analytically  a  value  for  the  steady  state  potential  drop  near  the  source  plate  for  monotonic 
potential  decay.  The  model  assumes  cold  ions,  injected  at  x  =  0  from  rest,  all  traversing  the  device, 
with  a  given  flux  density  T,,  independent  of  x.  The  warm  electron  density  is  assumed  to  be  governed 
by  a  Boltzmann  factor.  For  convenience,  the  source  is  at  zero  potential,  $(0)  =  0 

The  analysis  assumes  time  independence,  and  uses  conservation  of  energy  and  ion  flux  to  set  up 
Poisson’s  equation.  At  this  step,  one  obtains  what  appears  to  be  a  maximum  ion  flux  into  a  neutral 
plasma,  occuring  at  an  ion  velocity  equal  to  sound  velocity.  For  the  next  step,  integrate  Poisson  once, 
and  apply  the  boundary  condition  of  zero  field  at  the  plasma  (x  -*  «).  These  steps  relate  the  ion  flux 
and  the  source  sheath  potential  drop.  The  steps  are  similar  to  these  of  Bohm  (1949)  or  Chen  (1973)  in 
their  well  known  calculauons  of  the  collector  sheath  potential  drop. 

The  ions  are  assumed  start  from  rest  at  x  =  0.  Their  energies  are  given  by 


(1) 


\  Mu?  =  -«D(x)  >  0  .  u,  = 


The  ion  flux  density  is  given  by 


T,  =  n,(x)  u,(x) 


(2) 


which  is  constant  with  x.  All  ions  are  accelerated,  starting  from  rest  The  ion  density  is 


_  T, 


r, 


(3) 


which  blows  up  at  the  source,  where  the  ion  velocity  vanishes.  This  implies  a  very  copious  source, 
similar  to  the  infinite  density,  zero  velocity  source  of  Child’s  Law  (1911). 

The  electron  density  is  taken  to  be  governed  by  a  Boltzmann  factor 


n,(x)  =  n„  exp  £  e  <h(x  )/kT,  j 


(4) 


n.„  is  the  electron  density  at  the  source  plate,  including  both  injected  and  returned  electrons.  This 
assumption  is  on  shaky  ground,  as  even  at  the  source  plate,  where  the  distribution  is  a  cut-off  Gaussian. 

Using  these  densities,  Poisson's  equation  is 


d2<t>(x)  e 


dx2 


to 


to 


n0,  exp 


e  <t> 
kT. 
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-2e 
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(5) 


For  convenience,  let 


— <r  <t>  n 


t--r  = 


tJcT. 


(6) 


(7) 
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r,  1  _  n  o, 

<r\  ~  vtf 


notv,<2 


is  the  normalized  ion  flux. 


Hence,  Poisson’s  equation  becomes,  with  =  r)  (q) 


•/  r-  -*  -’I 

i\  =Fi\  ~e 


INSIGHT  INTO  ION  FLUX  and  ION  VELOCITY 


At  this  step,  without  employing  any  boundary  conditions  on  the  field  at  x  =  0  or  at  x  =  , 

we  can  find  that  value  of  potential  t)  and  flux  F  at  which  t| '  changes  sign,  that  is,  ti  '(T|,F)  =  0  at 


F(n'  =  0)  =  vn 


This  reLauon  Ls  plotted  in  Figure  2.  The  maximum  value  of  flux  is  found  at 


ti  =  0.5 


=  0.4288 


T,  =  =  0.6065  v, 

mu  yg 


V,  =  v2T|  V,  =  V, 


These  statements  say  that  the  maximum  ion  flux  into  a  neutral  plasma  occurs  with  a  potential  drop  from 


the  source  of  r|  =  0.5  (e  /  kT,  =  -0.5)  which  accelerates  the  ions  to  sound  velocity. 
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Figure  2.  Locus  of  normalized  ion  flux  density  F  as  a  function  of 
potential  11(4),  with  %  located  at  zero  net  charge  density, 
r|"(4)  =  0.  The  ion  velocity  v,  (£)/v,=  n'2t)~(4)  is  also  noted 


■Sv^’’CVL,vV-'r-.'C--Cs'' 


men 


•  -*-•  ,*»  *'•  b*'  b‘-  V  *.  .**  , 

,  •  -  -  «  ■  *  «  .  *  «  1 ■  ,  - 


FIRST  INTEGRATION;  USE  ZERO  FIELD  AT  PLASMA 


Integrate  Equation  (10)  once  by  multiplying  by  T|';  integrate  |  dt,  ,to  obtain 


-(TlO2 


4 

4 

=  2  F  T)1'4 

+  e 

0 

0 

(16) 


The  boundary  conditions  at  the  left  side  are  :  T)  =  0  at  4  =  0,  i\'  (0)  =  ti'0.  At  the  right  side  we  take 
4  to  be  in  the  plasma  region  where  we  postulate  zero  electric  field,  T]'  (4)  =  0.  Hence, 


4-oi'o)J  =  i  -  [«'n  +  2F t\,/4] 


(17) 


We  require  the  left  hand  side  to  be  zero  or  positive,  for  real  solutions;  hence, 

(18) 

where  ri  is  measured  at  sheath  edge,  4.  at  the  plasma  potential.  Tip. 


This  inequality  may  be  rearranged  in  order  to  be  read  as  a  limitation  on  the  ion  flux,  (for  tj'  (4) 


0). 


F  <  (l  -  e~'')  /  (2Vti)  =  G  (ti)  (19) 

For  small  T),  G(q)  =  /  2  ,  for  large  t),  G  Cn)  =  1/(2  vrj)  ;  in  between  these  values  G  (T|)  has  a  max¬ 

imum  value  of  Gm„  =  0.3191  which  occurs  at  t|  ■  r|c  =  1.25643 

Taking  F  <  0.3191  s  Fc 

T,  <0.3191  <2n0,  v,  =  0.45 1 3  n„  v,  =  (20) 

When  the  ion  flux  T,  <  I-*  (F  <  Fc),  then  there  are  two  values  of  t|  (t),  and  r|2  ,  see  Figure  3)  at 
which  F  =  GO)),  meaning  r|'(0).  In  between,  T),  <  t)  <  t|2,  F  <  G  (t))  so  that  (Tf2)  >  0.  Outside,  r|  < 
and  T|  >  t)2,  F  >  G  0D  so  that  (q'2)  <  0,  taken  as  an  indication  of  a  solution  that  does  not  fit  our 
assumptions  (possibly  a  me -dependent,  possibly  unstable,  possibly  multivalued  v,  ,  or  other). 


1.25643 


7)  0.5  1.0  T)c 


2.0 


Figure  3.  Flux  function  G(ri)  locating  regions  where  the  source 

field  squared  changes  sign.  There  is  a  maximum  ion  flux  value 
(Fc  =  0.3191)  of  a  critical  potential  (r{c  =  1.25643). 


For  ri,  <  r|c .  with  F  <  G  (q),  (q7(0))2  £  0,  Equation  (10)  gives  T|"i  (£)  <  0,  all  of  which  fits  the 
expected  q(£)  in  Figure  1.  However,  for  q2  >  qc ,  F  =  G  (tj)  q"2(£)  >  0,  which  does  not  fit  Figure  1. 

As  F  is  increased  toward  Fc ,  T],  and  q2  coalesce  to  qc  ,  where  F  =  Fc ,  q  =  qc ,  q70  =0. 

The  maximum  value  of  G(q)  is  found  at  G'(q)  =  o,  which  is  at 

1  +  2qc  = 

solved  by  qc  =  1.25643.  This  is  the  same  relation,  for  F  =  G  (q),  for  13"  =  0.  Hence  at  qc,  q"(ij)  =  0 
and  r|'(4)  =  0. 

When  an  F  >  Fe  (T,  >  T^)  is  imposed,  it  is  postulated  that  too  many  ions  are  injected  relative  to 
the  electron  injection,  that  the  potential  near  the  source  rises  and  reflects  some  ions  at  a  virtual  anode. 
This  is  an  easy  physical  concept,  but  is  not  allowed  here  analytically.  The  ions  here  are  cold,  all  with 
the  same  velocity;  the  analysis  needs  changing  for  virtual  anode  formation,  with  a  finite  ion  velocity 
distribution  which  allows  for  some  ion  flow  to  the  plasma  and  for  ions  returning  to  the  source. 

COMBINING  ZERO  NET  DENSITY  and  ZERO  FIELD  AT  PLASMA 

Let  us  combine  the  zero  net  charge  density  (13"  (£)  =  0)  and  the  zero  field  (13'  (£)  =  0)  informa¬ 
tion.  We  do  this  by  plotting  F  for  13"  (£)  =  0  (Figure  2)  and  F  for  Tj  (%)  =  0  and  q'  (0)  =  0  (Figure 
3)  on  the  same  axes,  as  appears  in  Figure  4.  Let  the  former  be  F  =  FD  =  e'1  Vq  and  the  latter  beF  = 
Fe  =  (1  -  e~n)  /  2  vij.  Then  let  us  ask  in  what  region  are  there  solutions  that  fit  the  model  at  Figure  1? 
With  the  areas  be  denoted  as  shown,  what  do  we  observe? 

areas  B  and  C  are  where  (q7  (0))2  <  0,  not  allowed  here. 

area  D  is  where  (q'  (0))2  >0  (all  right)  but  q"  (£)  >  0,  which  does  not  fit  Figure  1. 
area  A  is  where  (q'(0))2  >  0  and  q"  (£)  <  0,  which  does  fit  Figure  1 
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Figure  4.  Combination  of  the  ion  flux  FD  associated  with  zero  net 
charge  density  at  the  plasma  (t|"(^)  =  0)  of  Figure  2  with 
the  ion  flux  FE  associated  with  zero  electric  field  at 
the  source.  Only  solutions  in  area  A  fit  the  assumed  t|(£) 
shown  in  Figure  1.  All  of  the  above  assumes  =  0, 
zero  field  at  the  plasma. 
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SECTION  D:  PLASMA-WALL  PHYSICS,  THEORY  AND  SIMULATION 


A.  Collector  Sheath,  Presheath,  and  Source  Sheath  in  a  Collisionless,  Finite  Ion 
Temperature  Plasma 

Lou  Ann  Schwager 

B.  Potential  Drop  and  Transport  in  a  Bounded  Plasma  with  Ion  Reflection 
at  the  Collector 

Lou  Ann  Schwager 

C.  Potential  Drop  and  Transport  in  a  Bounded  Plasma  with  Secondary  Electron  Emission  at 
the  Collector 

Lou  Ann  Schwager 

Final  reports  with  the  above  titles  are  in  preparation  for  ERL  memos  and  journal  submissions. 
The  final  titles  may  be  altered  slightly. 

D.  Magnetized  Sheath;  Kelv  in-Helmholtz  Instability;  Long-lived  Vortices 
Dr.  Kim  Theilhaber 

A  movie  has  been  made  displaying  the  long-lived  vortices  in  detail.  An  ERL  memorandum  has 
been  prepared  to  accompany  the  movie. 
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SECTION  HI:  CODE  DEVELOPMENT  AND  SOFTWARE  DISTRIBUTION 


Codes  which  we  have  developed  or  modified  here  are,  generally  speaking,  available  for  the  ask¬ 
ing.  We  are  making  arrangements  with  our  Industrial  Liaison  Program  to  handle  distribution,  at  cost,  of 
listings,  tapes  and/or  diskettes.  In  a  later  report  we  will  provide  a  list  of  codes  and  related  material 
available  and  how  to  obtain  such.  If  you  use  those,  we  simply  ask  that  you  acknowledge  the  source 
and  also  send  us  your  results  as  published  in  reports  or  journals. 

(Not  all  of  the  codes  mentioned  in  our  reports  are  available  for  distribution;  some  a'e  single¬ 
purpose,  single-user,  too  stylized  for  outside  use.) 

A.  A  Relativistic  Monte  Carlo  Binary  Collision  Model  for  Use  in 
Plasma  Particle  Simulation  Codes 

R.  J.  Procassmi,  C.  K.  Birdsall,  E.  C.  Morse  and  B.  I.  Cohen 

This  report  is  to  be  issued  as  ERL  Memorandum  No.  UCB/ERL  M87-24,  May  14,  1987.  The 
introduction  follows. 

B.  Performance  and  Optimization  of  Direct  Implicit  Time  Integration 
Schemes  for  Use  in  Electrostatic  Particle  Simulation  Codes 

R.  J.  Procassmi,  C.  K.  Birdsall  and  E.C.  Morse 

and  B.  1.  Cohen  (Lawrence  Livermore  National  Laboratory) 
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A  Relativistic  Monte  Carlo  Binary  Collision 
Model  for  U  e  in  Plasma  Particle  Simulation  Codes 


R.  J.  Procassim,  C.  A.  Birdsall  and  E.  C.  Morse 

and 

B.  I.  Cohen  (Lawrence  Livermore  Sat  tonal  Laboratory) 

1.  Introduction 

Particle  simulations  of  plasma  physics  phenomena  employ  far  fewer  particles 
than  the  systems  which  are  being  simulated,  owing  to  the  limited  speed  and  memory 
capacity  of  even  the  most  powerful  supercomputers.  In  practice,  a  system  containing 
a>  many  as  I02CI  particles  is  studied  with  upwards  of  104  simulation  particles.  If  the 
simulation  consists  of  point  particles  in  a  gridless  domain,  then  the  combination  of 
the  small  number  of  particles  in  a  Debye  sphere  and  the  possibility  of  zero-impact- 
paranfter.  large-angle  scattering  results  in  a  significant  enhancement  of  fluctuation 
phenomena  such  as  collisions. 

In  contrast,  a  gridded  simulation  has  finite-size  particles,  due  to  the  weight¬ 
ing  or  interpolation  of  the  particles  onto  the  grid  and  vice  versa  for  the  forces. 
These  finite-size  particles  are  tenuous  in  nature,  being  able  to  pass  through  each 
other.  Therefore  the  interpart icle  force  does  '  ot  become  singular  as  the  distance 
of  separation  of  the  particle  midpoints  goes  to  zero,  as  it  does  for  point  particles.1 
Tin’s  means  that  (he  average  scattering  angle  for  finite-size  particle  interactions  is 
much  smaller  than  that  for  point  particle  interactions,  hence  the  collisionality  for 
finite-size  particles  is  reduced  relative  to  that  for  point  particles. 

Accommodating  collisional  processes  in  a  simulation  may  be  difficult  because 
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of  disparate  time  scales.  A  comparison  of  the  relevant  physical  time  scales  of  the 
system  that  is  being  simulated  usually  yields  a  large  range  of  values.  For  instance, 
the  grid-cell  transit  time  is  usually  several  orders  of  magnitude  smaller  than  the  90 
scattering  time.  Much  of  the  physical  phenomena  of  interest  in  the  simulation  are 
due  to  these  long-time-scale  collisional  processes2,  but  short-time-scale  processes 
(such  as  particle  bounce  times  in  a  mirror  or  tokamak)  must  be  adequately  resolved 
if  the  plasma  dielectric  response  and  the  plasma  potential  are  to  be  accurately  deter¬ 
mined.  This  leads  to  important  constraints  on  the  modeling  of  binary  collisions  in 
our  electrostatic  particle  simulation  code  TESS  (Tandem  Experiment  Simulation 
Studies). 

In  the  following  section  we  outline  the  physics  and  operation  of  the  binary 
collision  model  within  the  electrostatic  particle  code.  Section  3  presents  the  results 
of  computer  simulations  of  velocity  space  transport  which  were  run  to  test  the 
accuracy  of  the  model.  Finally.  Section  4  discusses  the  timing  statistics  for  the 
collision  package  relative  to  the  other  major  physics  packages  in  the  code,  as  well 
as  recommendations  on  the  frequency  of  use  of  the  collision  package  within  the 
simulation  sequence. 


B.  Performance  and  Optimization  of  Direct  Implicit 
Time  Integration  Schemes  for  Use  in  Electrostatic 
Particle  Simulation  Codes 

R.  J.  Pnocassmi.  C.  A.  Bird  soil  and  E.  C.  Morse 

and 

B.  I.  Cohen  (Lawrence  Livermore  Xational  Laboratory ) 

1.  Introduction 

Implicit  time  integration  schemes  allow  for  the  use  of  larger  time  steps  than 
conventional  explicit  methods,  thereby  extending  the  applicability  of  kinetic  parti¬ 
cle  simulation  methods.  This  paper  will  describe  a  study  of  the  performance  and 
optimization  of  two  such  direct  implicit  schemes,  which  are  used  to  follow  the  tra¬ 
jectories  of  charged  particles  in  an  electrostatic,  particle-in-cell  plasma  simulation 
code.  The  two  schemes  that  are  studied  were  developed  and  analyzed  bv  Langdon1. 
and  by  Cohen.  Langdon  and  Friedman2.  The  direct  implicit  method  that  was  used 
for  this  study  is  an  alternative  to  the  moment -equation  implicit  method  developed 
by  Mason3  and  Denavit4. 

The  organization  of  this  paper  is  as  follows:  Section  2  will  present  the  formu¬ 
lation  of  the  two  schemes,  as  well  as  the  incorporat ion  of  these  schemes  into  the 
particle  simulation  code.  Section  3  will  present  the  results  of  a  numerical  stability 
study  of  these  schemes.  This  study  is  based  upon  the  energy  conservation,  or  lack 
thereof,  of  a  freely  expanding  plasma  slab  in  various  regions  of  u.'peAt  .  Ax  Xor 
parameter  space,  where  -Jp,.  is  the  electron  plasma  frequency.  is  the  electron 
Debye  length.  At  is  the  time  step  and  Ax  is  the  grid  spacing.  Finally.  Section  4 
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summarizes  the  results  of  this  study,  and  gives  suggestions  for  application  of  the 
results  to  the  choice  of  a  stable  operating  point  for  other  types  of  plasma  systems. 

2.  The  Direct  Implicit  Electrostatic  Particle  Code 

2.1  Direct  Implicit  Time  Integration 

The  two  implicit  time  integration  schemes  that  were  used  for  this  study  were 
chosen  for  their  desirable  properties.  The  most  important  of  these  are  ij  the  re¬ 
laxation  of  -.'A t  constraints  on  the  stability  of  the  method.  11)  strong  damping  of 
high-frequency  modes  for  which  *A t  ^  1.  and  in:  the  second-order  accuracy  of  the 
methods  in  the  simulation  of  low-frequency  phenomena  for  which  ^'Ai1  1.  These 

properties  were  originally  studied  through  the  normal-mode  analysis  of  an  elec¬ 
trostatic.  cold-plasma  oscillation  in  which  the  electrons  undergo  simple-harmonic 
oscillations.2  The  present  study  deals  with  one  scheme  from  each  of  the  C’  and  D 
classes  of  schemes  that  are  discussed  in  2  .  We  will  now  briefly  review  the  descrip¬ 
tion  of  the  C’]  and  D j  schemes. 

2.1.1  The  C  1  Scheme 

The  first  multistep  scheme  that  shall  be  considered  for  use  in  the  time  integra¬ 
tion  of  the  particle  trajectories  is  of  the  form 

xn  -  iC'oa,,  -  Cha.,-1  -  ...  -  0-2a„-k-2  )Ar  -  x„  ( la  1 

a»-i  =  (x|,  -  2x^_!  -  x|,_2)  A/2.  tVn 

where  x  is  the  particle  displacement,  x'  is  the  leapfroe  particle  displacement  and  a 
is  the  acceleration.  The  subscripts  indicate  the  time  step  level  of  the  variable,  li  x' 


is  eliminated  from  the  equations  we  get 


(x„  -  2xn _ i  ~  xn_2  )  At2  —  a„_ i  -  C0(a„  -  2an_:  -  a„_2)"- 

f'i  *  an-  1  -  2a„_2  -  a„-s)  - 

(  vlan-k  ~  -an  -  k  -  1  ~  an-k-2)  —  •  ■  (2) 

This  method  is  second-order  accurate  for  small  tune  step.  The  so-called  C\  scheme 
has  the  coefficients  C,  =  0  for  s  _  2.  A  stability  analysis  for  this  scheme  was 
performed  by  varying  C\>  and  C\  independently,  in  an  effort  to  determine  the  min¬ 
imum  value  of  c  =  exp\  -  )u.Af  i  for  the  least  damped  mode  of  a  simple-harmonic 
oscillation  at  a  frequency  -cr>.  The  values  oi  C' o  and  (\  at  which  these  minima  occur 
are  C o  -  0.302  and  C\  -  0.040.  corresponding  to  an  optimized  scheme  for  which 
z  --  0.5  as  o.'P2!  —  tx.  Therefore,  the  "optimized  C j  scheme’’  is  given  by 

(xn  -  2x„_i  -  xn_2l-lr  =an_!  *  0.302C an  -  2a„_i  -  a„-2)- 

0.0401  an_x  -  2a„_2  -  an  3).  (3) 


2.1.2  The  Di  Scheme 

The  next  multistep  scheme  that  was  u>ed  was  derived  from  a  classic  second- 
order  stiffly  stable  scheme  for  the  solution  of  a  first -order  ordinary  differential  equa¬ 
tion  of  the  form  dy  dt  —  f'y.t)-  This  multistep  method  is  of  the  form 

y,  ~  i  --  1  4  3  \yK  ~  <  1  3  \y.,  -  i  =  ( 2  3  )/r. _  j  Af .  ( 4  i 

The  form  of  tliis  equation  motivates  consideration  of  the  more  general  class  of 
second-order  -rhemes  which  u-e  dx  dt  --f  v  and  </v  df  =  a  to  produce 


In  order  to  make  this  method  second-order  accurate,  it  was  found  that  the  coeffi¬ 
cients  .  ii i .  d j . .  .  .  have  to  be  chosen  to  satisfy  the  following  two  criteria 

1  =  —  d  i  —  d;  —  .  .  .  I  (ia  ) 

and 

0  =  d<>  —  2d  i  —  3d 2  —  . .  •  .  i  (ib  ! 

The  simplest  of  the  class  D  algorithms  has  d,  —  0  for  s  -f  2.  such  that  from  |6ai 
and  (bh)  we  find  that  dr,  =  2  and  dj  =  —  1.  This  "D\  scheme"  is  then  given  by 

-  2x„_1  -  oxIt  -  lx„_x-xn_2-  ' 

2.1.3  General  Results  and  Properties 

It  was  sho”  .  in  2  that  relaxation  of  the  wpAf  stability  constraint  in  plasma 
simulation  required  lie  use  of  an  implicit  time  integration  scheme.  This  means  that 
temporal  resolution  of  the  plasma  frequency  is  forfeited  to  simulate  lower  frequency 
phenomena  efficiently.  This  is  particularly  useful  if  the  phenomenon  of  interest 
for  the  simulation  is  on  a  much  slower  time  scale  than  t hat  of  plasma  oscillations. 
Therefore,  a  tim<-  step  that  encompasses  several  plasma  oscillation  periods  may  be 
used  in  lone  time  scale  simulations  while  the  scheme  remains  numerically  stable. 
The  coefficient  >  that  are  med  in  the  time  integration  schemes  determine  the  decree 
of  implicitness,  ’he  degree  of  high-frequency  mode  dissipation,  and  the  stability 
of  the  scheme.  When  1.  the  high-frequency  oscillations  are  distorted, 

even  though  the  scheme  remains  stable.  Therefore,  the  scheme  should  damp  out 


these  distorted  inodes,  while  not  affecting  the  low-frequency  modes  which  are  of 
interest.  The  c  for  the  least  damped  simple-harmonic  oscillator  normal  mode  U 
plotted  as  a  function  of  *02A t~ ,  for  both  the  Cfi  and  D\  schemes  in  Figure  1.  The 
figure  shows  that  the  damping  for  the  D\  scheme  is  stronger,  and  begins  at  smaller 
values  of  At  ,  than  that  for  the  C'j  scheme.  This  is  due  to  the  fact  that  the  only 
acceleration  used  in  (7)  is  a„._x-  indicating  that  the  D ]  scheme  attempts  to  force 
charge  neutrality  in  only  one  time  step  as  A/  -  x.  In  contrast,  the  C'j  scheme 
uses  acceleration  data  from  several  previous  time  steps.  Therefore,  the  Cj  scheme 
requires  that  more  previous  time  step  data  be  stored,  than  does  the  Dx  scheme.  The 
storage  of  this  data  could  become  an  important  issue  as  the  number  of  particles  used 
in  the  simulation  grows.  The  accuracy  of  the  two  schemes  at  low-frequencies  was 
discussed  in  2  .  The  weak  damping  of  low-frequency  oscillations  indicated  in  Figure 
1  yields  high  accuracy  simulations  in  the  "explicit"  regime  o.'0A/  <  1. 

2.2  Electrostatic  Field  Equations 

In  addition  to  the  implicit  time  integration  of  the  particle  trajectories,  the 
calculation  of  the  field  quantities  represents  the  other  major  implementation  issue. 
The  new  position  xn-i  was  shown  to  depend  upon  the  new  acceleration  an.j  in 
the  implicit  time  integration  schemes.  In  an  implicit  electrostatic  particle  code, 
this  acceleration  is  due  to  the  new  electric  field  E,,.!.  However,  this  field  is  not 
yet  known,  since  it  depends  upon  the  charge  density  p„  _ !  and  the  particle  position 
xn^i-  Therefore,  we  have  a  large  system  of  coupled,  nonlinear  particle  and  field 
equations  that  must  be  solved. 
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Figure  1.  The  absolute  value  of  z  =  exp(-t-zAf)  of  the 
least  damped  simple-harmonic  oscillator  normal  mode  versus 
in  the  C\  scheme  (Eq.  (3):  dashed  line)  and  the  D\ 
scheme  (Eq.  (7):  dashed-dotted  line). 
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One  early  method  that  was  used  to  calculate  the  field  quantities  at  the  advanced 
time  level  involved  the  solution  of  a  set  of  coupled  field  and  fluid  equations3  ,4. 
This  is  known  as  the  implicit  moment -equation  method.  In  contrast,  the  direct 
implicit  method  uses  a  linearization  of  the  particle-field  equations  to  calculate  the 
new  electric  field  En_i. 

In  general,  the  implicit  time  integration  schemes  described  above  give  the  new 
particle  position  as 

xn_i  =  iAt2an_1  4-  x,  (8) 

where  0  <  3  1  is  a  constant  of  the  time  integration  scheme  and  x  is  the  explicitly 

computable  "free-streaming''  particle  position  .  which  is  independent  of  an_f  In 
an  electrostatic  code5,  an_!  is  due  only  to  En_x. 

The  new  position  xn-i  may  be  thought  of  as  the  explicitly  computable  position 
x  plus  a  displacement  6x  =  3 At2 a„_i.  The  charge  density  is  then  taken  as  the 
sum  of  p  depending  on  the  positions  of  X.  plus  a  linearized  correction  term  due  to 
the  displacing  of  the  particles  by  the  amount  6x  —  xn-i  -  x,  which  is  equal  to 

Sp  ~  -  V  •  [p(x)6x(x)'.  (9) 

The  linearized  displacement  term 

£x(x)  =  JAt2gE„_1  m  (10! 

then  eives 

£p(x)  =  -V  •  'y(x)E„*i  ,  (11 ) 
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where  the  effective  susceptibility  is  defined  bv 

\(x)  =  3  qp(x)  m  =  3^'2Al2.  (12) 

(These  expressions  have  been  generalized  to  the  relat ivisitic  case  in  our  computer 
code  TESS.  which  stands  for  Tandem  Experiment  Simulation  Studies).  Poisson's 
equation  is  modified  by  the  addition  of  t His  susceptibility  to  become,  in  rationalized 
cgs  units 

V.E„^  V.(.xEn+1)  (13) 

or 

-T  ■  T  ~  \  Vc>n-1  =  P-  (14) 

This  form  of  the  field  equation  depends  only  upon  the  particle  position  xn-i.  as  op¬ 
posed  to  the  implicit  moment -equation  method  which  also  requires  particle  velocity 
information5 . 

The  TESS  code  uses  a  simplified  finite  difference  form  of  (  14).  which  in  the 
absence  of  local  spatial  smoothing  is  given  by:5 

-”1  ~  ~  .2  -  (\j  ~  \j-iY.Oj  -  T  -  Xj-T.Oj-j  =  Pj  Ai2  (15) 

where 

\j  =  i\j  ~  (16) 

is  the  midpoint  effective  susceptibility. 

In  the  simulation  of  long  time  scale  phenomena,  it  is  desirable  to  have  a  value 
of  •  1.  which  means  that  1  —  \  1,  such  that  the  field  equation  is 
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substantially  different  than  Poisson's  equation,  which  is  the  explicit  limit  of  (1-4)  for 


3.  Direct  Implicit  Method  Performance  Study 

The  performance  of  the  Cj  and  D\  time  integration  schemes,  as  a  function 
of  location  in  .  A.r  parameter  space,  was  analyzed  for  the  collisionless 

expansion  of  a  one-dimensional,  sharp  boundary  plasma  slab  into  a  vacuum.  The 
free  expansion  of  such  plasma  slabs  has  been  previously  discussed  by  Denavit4,6. 
including  a  comparison  of  particle  simulation  results  to  a  self-similar  solution  which 
assumes  charge  neutrality,  isothermal  electrons  and  cold  ions.  As  the  plasma  ex¬ 
pands  into  the  vacuum  due  to  the  velocity  of  the  particles,  a  rarefaction  wave 
propagates  into  the  plasma  slab  at  the  ion-acoustic  speed  e3.  The  plasma  density 
behind  the  wavefront  falls  off  exponentially,  and  the  ions  are  accelerated  outward 
according  to  the  velocity  profile 

Vf ,  =  c,-  (x  -  x0)/t  (17) 

where  vj,  is  the  velocity  of  the  expanding  ion  front,  xq  is  the  location  of  the  initial 
density  discontinuity  and  x  is  the  location  of  the  ion  front  at  time  t.  Figure  2  shows 
the  electron  and  ion  phase  space  plots  for  a  typical  expanding  slab.  The  dashed  line 
on  the  ion  phase  space  plot  (Figure  2b)  represents  the  self-similar  ion  front  velocity 
as  given  by  (17).  The  agreement  between  the  self-similar  and  particle  simulation 
ion  front  velocities  shown  in  Figure  2  was  seen  earlier  in  6  . 

Freely  expanding  plasma  slab  simulations  were  made  for  various  combinations 
of  the  important  simulation  parameters  and  &x/\£)e  in  the  range 
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0.1  _  u.'pt  At  ^  200.0,  0.1  <  Aj  Ape  <  100.0.  excluding  the  region  for  which 
r,/;,  Af  -A.r  •  2.  Several  input  parameters  were  the  same  for  each  of  the  simu¬ 
lations.  These  include  the  number  of  particles  (2048  of  each  species),  the  mass 
ratio  (no  m,  =  900u  the  temperature  ratio  ( Te  T,  ~  10).  the  system  half-width 
l  L  =  100cm  n  the  initial  plasma  slab  half-width  [L,iab  —  30cm  I  and  the  mesh  size 
!  As  —  L  12s  -  0.?M3cmi.  The  values  of  and  A.r  were  the  same  for  each  of 
the  simulations,  while  the  temperature  of  the  species  (and  hence  Ap)  and  At  were 
varied. 

In  an  attempt  to  quantify  the  performance  of  the  two  time  integration  schemes 
in  various  regions  of  At  .  A.r  \[,.  parameter  space,  we  employ  the  quantity 
A  E  E,  i  A.  where  AF  is  the  chan  tie  in  the  total  plasma  energy  (kinetic  plus 
electrostatic  !  through  .V  time  steps,  normalized  to  the  initial  value  of  total  plasma 
energy  Figure  3  shows  a  typical  time  history  of  the  total  plasma  energy  for  a 
simulation  that  resulted  in  numerical  heating  of  the  plasma.  The  freely  expanding 
slab  should  exactly  conserve  energy,  since  the  kinetic  energy  of  the  hot  electrons  is 
simply  transferred  to  the  cold  ions  through  the  self-electric  field,  resulting  in  the 
acceleration  of  the  ions  which  was  observed  in  Figure  2b.  However,  since  TESS  uses 
finite  difference  equations  and  the  implicit  algorithms  described  here  are  inherently 
dissipative,  energy  is  r-nt  necessarily  conserved.  Some  possible  causes  or  mechanisms 
of  non-conservation  of  energy  include  >  •  stochastic  self-healing  effect-  r*.  u1  self¬ 
heat  i Tig  effect  s  arising  from  the  interact  i<  >n  <  •!  -pat  ial  aliases,  which  are  a  consequence 
of  the  discrete  nature  oi  the  grid,  with  the  velocity  distribution  function  <>i  the 
particles10  and  in)  numerical  heat  mg  and  cooling  due  to  the  dissipative  nature  of 
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the  implicit  time  integration  scheme2 .  This  section  presents  the  results  of  a  study 
of  t he*  dependence  of  energy  conservation  on  -c?t.A t  and  Ax  Xp^.  using  the  C\  and 
D i  implicit  time  integration  schemes. 

A  set  <>t  3'>  simulations  were  run  for  various  combinations  of  A/  .  A.r  \/,( 
which  span  the  aforementioned  space.  These  runs  were  made  using  the  Ch  and  /)j 
schemes.  both  with  and  without  one-pass,  local,  self-consistent  digital  smoothing 
of  (14*  as  described  in  detail  in  Section  3.3  of  5  .  The  spatial  smoothing  that  was 
applied  was  of  a  simple,  three-point.  (1.2.11  form.  It  was  therefore  also  possible  to 
study  the  effect  of  spatial  smoothing  on  the  non-conservation  of  energy.  The  energy 
convert  at  ion  results  for  these  runs  are  given  in  Tables  1  and  2  for  the  (\  scheme 
with  and  without  spatial  smoothing  and  in  Tables  3  and  4  for  the  I>i  scheme  with 
and  without  smoothing  respectively.  The  data  from  these  tables  is  shown  in  the 
contour  plots  of  Figures  4  through  7. 

The  zero  contour  line  on  these  figures  indicates  exact  energy  conservation,  while 
positive  i  negative i  contour  lines  indicate  numerical  heating  ( cooling  i  of  the  plasma, 
i  The  positive  contour  lines  are  solid,  the  negative  contour  lines  are  dashed,  the  zero 
contour  line  i-  the  bold  solid  curve  and  the  c.»(. A t  Ax  —  1  line  is  the  bold  dashed 
line i .  I  he  filied  circles  show  the  locations  in  .c;,(At  .  A.r  \p,  parameter  space  at 
which  *  iv  simulations  were  run.  The  coarse  nature  of  the  contour  lines  is  due  to 
the  limio-d  number  ot  runs  in  the  simulation  set.  The  space  could  be  extended  to 
larger  value-  o)  the  two  parameters,  but  cost  of  each  run  increases  rapidly  as 
Ax  A i ,,  is  increased,  due  to  the  decreasing  particle  thermal  velocity,  which  in  turn 
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Table  1 


Energy  Conservation  for  the  C\  Scheme  without 
Spatial  Smoothing  in  xpt  At  ,  Aj  \De  Parameter  Space 
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Table  1  (continued) 
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Figure  4.  Contour  plots  of  (A E/E0)/N  as  a  function  of  u>p,A t 
and  At/A p,  for  the  C,  Scheme  without  Spatial  Smoothing. 
The  zero  contour  indicates  exact  conservation  of  energy. 


Table  2 


Energy  Conservation  for  the  C' \ 
Spatial  Smoothing  in  *\„  At  ,  A.r  \[,t 

Scheme  with 
Parameter  Space 

A/ 

A.r  \lu 
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Table  2  (continued) 
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A./'  \[>. 
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re  5.  Contour  plots  of  (A E/E0)/N  as  a  function  of  u^A 
Ax/\[)e  for  the  C\  Scheme  with  Spatial  Smoothing.  The 
contour  indicates  exact  conservation  of  energy. 


requires  a  greater  number  of  time  steps  for  a  uniform  expansion  of  the  slab  among 
all  of  the  runs. 

Explicit  particle  simulation  codes  are  usually  limited  to  operation  in  the  lower 
left  hand  region  above  the  bold  dashed  line.  The  direct  implicit  method  greatly 
extends  this  region  of  operation.  The  figures  indicate  that  it  is  possible  to  use 
large  values  of  At  and  A.r  \pf  which  result  in  negligible  or  no  loss  of  energy 
conservation.  Therefore,  one  may  increase  the  cost  effect iveness  of  the  simulation 
by  employing  a  time  step  which  is  larger  than  the  electron  plasma  period,  while 
accurately  simulating  the  low-frequency  phenomena  of  interest,  such  as  MHD  (in 
implicit  electromagnetic  codes :  and  particle  transport.  It  must  be  stressed  that  the 
results  presented  here  are  applicable  for  either  periodic  plasma  systems  or  systems 
in  which  t lie  plasma  is  not  in  contact  with  a  wall.  The  effect  of  potential  sheaths 
near  a  surface  may  yield  different  conclusions  regarding  regions  of  numerical  heating 
and  cooling  than  those  presented  here. 

The  contour  plots  for  the  C\  scheme.  Figures  4  and  5.  show  that  energy  is 
conserved  up  to  maximum  values  of  A.r  \p(  =:  3  and  ^peAt  1.  For  larger 
values  of  these  parameters,  there  is  numerical  heating,  while  smaller  values  lead  to 
numerical  cooling  of  the  plasma.  The  addition  of  local  spatial  smoothing  allows 
the  zero  contour  line  to  extend  to  slightly  larger  values  of  the  two  parameters. 
Smoothing  also  reduces  the  numerical  heating  at  large  values  of  A.r  \^,r  by  about 
a  factor  oj  4.  The  data  points  at  u.\,tAt  —  200,0  .  At  \p,  =  100.0  are  subject  to 
very  rapid  numerical  heating  with  the  C' j  scheme.  Spatial  smoothing  increases  the 


heating  rate  at  this  point,  which  is  the  opposite  effect  to  that  seen  in  other  regions 
of  the  parameter  space. 

In  contrast  to  the  contour  plots  for  the  Cj  scheme.  Figures  6  and  7  indicate 
that  there  are  regions  in  the  space  with  large  values  of  -xpeA t  and  Ax,  \pe  for 
which  the  energy  of  the  system  is  conserved  well  when  the  Dj  scheme  is  used.  The 
zero  contour  line  lies  approximately  parallel  to  the  r(fceA t  Ax  =  1  (bold  dashed) 
line,  and  extends  out  to  wpeA t  =  30  .  Ax  \oe  —  100  following  a  path  described 
approximately  by  r,/,eAf/Ax  0.30  ~  0.1 3  for  uvpeA t  >  1.  This  behavior  may 
continue  to  larger  values  of  the  two  parameters,  which  would  allow  for  very  large 
time  step  simulations  in  which  numerical  heating  mid  cooling  are  not  significant. 
Note  that  there  is  a  very  rapid  rate  of  change  in  i  A E  Ec,  I  A  about  the  zero  contour 
line,  especially  for  A t  >  1.  Spatial  smoothing  yields  approximately  a  factor  of  3 
reduction  in  the  numerical  heating  at  large  values  of  Ax.  Aoe.  The  smoothing  also 
results  in  numerical  heating  for  small  values  of  the  two  parameters,  which  is  the 
opposite  of  the  results  obtained  without  any  spatial  smoothing. 

4.  Conclusions 

The  C' i  and  D\  direct  implicit  time  integration  schemes  have  been  shown  to 
extend  the  region  of  operation  for  which  the  toted  plasma  energy  is  conserved  beyond 
that  which  is  applicable  for  an  explicit  time  integration  scheme.  The  C j  scheme 
wax  found  to  have  a  limited  region  of  parameter  space  over  which  energy  is  well 
conserved.  The  D j  scheme  possesses  a  much  larger  region  of  parameter  space  for 
good  energy  conserving  operation.  In  fact,  for  the  range  of  parameter  space  0.1  £ 
-xp, At  _  200.0.  0.1  y  Ai  Apf  _  100,0.  the  D\  scheme  was  shown  to  conserve 


Table  3  (continued) 
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Scheme 


Contour  plots  of  (AE/  E0)/N  as  a  function  of  u ip, 
ip,  for  the  D\  Scheme  without  Spatial  Smoothing 
contour  indicates  exact  conservation  of  energy. 


Table  4 


Energy  Conservation  for  the  D}  Scheme  with 
Spatial  Smoothing  in  ,  Aj.'A D,  Parameter  Space 
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energy  well  along  a  path  described  by  rf),(A<  /Ai  0.30  ±0.13.  Numerical  heating 
was  found  for  large  values  of  A/  \pe  with  each  of  the  schemes.  The  amount  of 
heating  in  a  given  region  of  the  space  was  found  to  be  comparable  between  the  two 
schemes,  with  the  Dj  scheme  heating  slightly  less  than  the  Cj  scheme.  Numerical 
cooling  of  the  plasma  system  was  found  for  the  lowest  values  of  the  parameters  in  the 
space,  except  when  the  D\  scheme  was  used  with  local  spatial  smoothing.  Spatial 
smoothing  generally  yielded  a  factor  of  3  to  4  reduction  in  the  amount  of  numerical 
heating,  and  a  slight  decrease  in  the  amount  of  numerical  cooling  relative  to  the 
same  simulation  run  without  smoothing.  The  smoothing  did  not  substantially  alter 
the  gross  topology  of  the  energy  conservation  plots. 

If  low-frequency  phenomena  are  of  interest,  as  opposed  to  high-frequency  ef¬ 
fects,  the  large  degree  of  high-frequency  damping  provided  by  the  D\  scheme  (see 
Figure  1)  will  allow  the  use  of  time  steps  that  encompass  several  electron  plasma 
periods.  The  *'pf  At  time  step  constraint  is  thereby  relaxed,  resulting  in  an  increased 
cost  effectiveness  of  the  simulation.  For  such  simulations,  the  D\  scheme  is  more 
desirable  for  use  than  the  Cj  scheme.  Since  the  energy  conservation  contour  plot 
for  the  D\  scheme  has  a  steep  gradient  about  the  energy  conserving  contour,  care 
must  be  taken  to  ensure  that  the  operation  point  is  chosen  as  close  to  the  energy 
conserving  contour  as  possible. 

The  results  presented  here  are  for  a  freely  expanding,  bounded  plasma  slab. 
The  conclusions  that  are  drawn  from  these  results  are  probably  applicable  to  other 
plasma  systems.  A  similar  performance  study  should  be  undertaken  for  systems 


in  which  the  plasma  is  in  contact  with  a  wall,  in  order  to  determine  the  numerical 
performance  of  such  simulations  using  these  direct  implicit  schemes. 
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Pierce  Diode  Instability  Simulations  with  External 
Inductance  and  Resistance  --  T.L.  Crystal  and  S.  Kuhn, 
V. C. Berkeley ,  Electronics  Research  Laboratory.  The  classical 
Pierce  diode  is  a  bounded  plasma  with  an  external  short 
circuit;  the  linear  instability  analysis1  of  the  classical 
equilibrium  was  verified  using  lDlv  particle  simulations  , 
Linear  analysis  of  the  extended  Pierce  diode  [meaning  the 
same  uniform  equilibrium  but  with  finite  external  circuit  R 
or  L ]  predicts  significant  departures  from  the  classical  case 
in  terms  of  instability  thresholds,  growth  rates,  and  the 
physical  nature  of  the  instability.  The  peculiarities  of  these 
extended  Pierce  diode  eigenfrequencies  and  corresponding 
field  profiles  have  now  also  been  examined  with  particle 
simulations. 

•  Work  supported  by  ONR  contract  N00014-85-K-0S09,  and 
by  Austrian  Research  Funds. 

1.  S.  Kuhn,  Fhys.  Fluids  27,  1B34  (1935). 

2.  T.L  Crystal  and  S.  Kuhn,  Fhys.  Fluids  28,  2116  (1935). 


Transport  across  the  Plasma-Wall  Sheath 
with  Secondary  Electron  Emission  and  Ion  Re¬ 
flection  .  L.A.Schwager  and  C  .K . Birdsall ,  U  .C . 
Berkeley*  --  Particle  and  energy  transport  from 
a  warm  plasma  to  a  floating  wall  is  obtained 
with  electrostatic  particle  simulations  in  Id. 
The  wall  may  emit  secondary  electrons  or  re¬ 
flect  ions.  For  both  boundary  effects,  we  find 
that  the  collected  electrons  carry  an  average 
energy  to  the  plate  equal  to  the  bulk  plasma 
temperature  (kT  per  electron)  even  though  Te 
drops  by  two  in  the  sheath.  As  the  ion  reflec¬ 
tion  fraction  increases,  the  reflected  current 
increases  the  plasma  density  which  then  in¬ 
creases  the  ion  energy  flux  to  the  collector. 
Even  with  an  ion-ion  two-streaming  instability, 
the  time-averaged  values  of  transport  results 
agree  with  time-independent  analysis.  Increas¬ 
ing  secondary  emission  beyond  where  En  reverses 
increases  potential  fluctuations  and  causes  de¬ 
parture  from  time-independent  theory,  ^  es¬ 
pecially  for  cool  secondaries. 

work  supported  by  DOE  contract  1-442427-22401. 
iG.D. Hobbs  and  J  .A .Wesson , Plasma  Physics , 9, 196 7 


Electrostatic  Effects  on  Confinement 
Outside  The  Separatrix  of  Field  Reverse 
Conf iQurati ons ,  S.  E.  Parker,  C.  K.  Birdsall,  and 
K.  Theilhaber,  UC-Berkelev* — A  two  dimensional 
particle  code  is  being  used  to  study  the  formation 
of  electrostatic  potentials  due  to  magnetic  peaks 
and  wells  with  an  injected  flux  of  well 
magnetized  electrons  and  weakly  magnetized  ions. 
Electrostatic  effects  may  explain  the  enhanced 
plasma  confinement  observed  in  the  open  field 
line  region  of  the  FRC1 .  Exper i menta 1  results 
have  shown  the  formation  of  large  potentials  due 
to  nonuniform  magnetic  fields.  This  type  of 

phenomena  is  being  modelled.  We  are  seeking  to 
show  that  magnetic  wells  cause  the  formation  of 
potential  peaks  and  magnetic  peaks  form  potential 
wells  when  ions  are  weakly  magnetized  compared  to 
the  electrons.  These  results  will  be  used  to 
interpret  the  electrostatic  potential  contour  of 
the  open  field  line  region  in  the  FRC. 

1.  L.  E.  Bteinhauer i  "Enhanced  Confinement 

Outside  the  Separatrix  of  PRC’s”,  Spectra 
Technology,  Inc.,  (1986). 

*  Work  supported  by  DOE  contract  FGO 3 -8 6 ER5 3 2 2 0 


Direct  Implicit  Particle  Simulation  of  Tandem 
Mirrors.*  R.  J.  PROCASSINI,  C.  K.  BIRDSALL, t  and 
B.  I.  COHEN,  Lawrence  Livermore  National  Laboratory — We 
developed  TESS,  a  one-dimensional  (axial),  relativistic, 
electrostatic  particle  code  to  simulate  thermal  barrier 
formation  and  ambipolar  potential  particle  confinement 
in  one-half  of  a  symmetric  tandem  mirror  plasma.  Our 
recent  advances  in  TESS  development  include: 
optimization  studies  of  two  direct  implicit,  finite 
difference  schemes  yielding  operating  regimes  of 
relative  numerical  stability  (i.e.,  no  numerical  heating 
or  cooling)  in  u  gAt  vs  Ax/l^g  parameter  space  for  each 
scheme;  calculations  (in  the  Ibsence  of  auxiliary 
heating)  of  ambipolar  potential  of  a  collisionless 
plasma  in  a  non-uniform,  tandem  mirror  magnetic  field 
with  non-reflecting,  no-axial-loss  boundary  conditions; 
incorporation  into  TESS  of  a  relativistic,  energy  and 
momentum  conserving,  Monte  Carlo  binary  particle 
collision  algorithm.  We  will  also  review  our 
implementation  of  algorithms  for  rf  and  neutral  beam 
auxiliary  heating,  axial  and  radial  particle  losses,  and 
plasma  fueling. 

♦Work  supported  by  the  USDGE  under  Cor.t.  W-7k05-ENG-48 . 
TUniversity  of  California,  Berkeley,  CA. 
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Structure  of  the  Crosse J-Fiold  Sheath  -  K.  Theil'naber 
and  C.K.  Birdsail,  University  of  California,  Berkeley. 

We  consider  the  equilibrium  and  time-dependent  beha¬ 
vior  of  the  plasma-wall  sheath  which  is  formed  when  the 
external  magnetic  field  is  parallel  to  the  wall. The  pre¬ 
sence  of  the  material  boundary  initially  creates  a  layer 
of  charge  depletion  near  the  surface,  of  thickness  compa¬ 
rable  to  the  ion  gyro-radius.  This  results  in  a  strongly 
nonuniform  electric  field  which  induces  a  sheared  EX& 
motion  of  ions  and  electrons  parallel  to  the  wall.  Physi¬ 
cal  arguments  as  well  as  analytical  and  numerical  studies 
of  a  similar  but  unbounded  configuration 1)  suggest  that 
this  flow  is  subject  to  the  Kelvin-Helmhol tz  instability, 
with  a  long-term  evolution  into  a  set  of  vortices  separa¬ 
ting  the  boundary  from  the  plasma.  We  shall  present  an 
analytical  model  for  the  instability,  and  the  results  of 
two-dimensional,  electrostatic  particle  simulations  of 
this  configuration,  with  emphasis  on  time-average  trans¬ 
port. —  DOE  Contract  FG03-S6ER33  220,  ONR  Contract  NO 00 14- 
85-K-0809 . 

(1)  P  .L.Pri tchett ,  F.V .Coroniti ,  J .Geophys .Res . , Vol .89 , 
No.Al,  p.168  (  198.)  . 


Bounded-Plasma  Physics  Workshop 
University  of  California,  Berkeley 
September  22-23,  1986 

Schedule  of  BPPW  Presentations 
Monday  22  September  1986 

8:30  Registration  and  parking  arrangements. 

Coffee,  tea,  orange  juice  and  rolls. 

9:00  C.K.  Birdsall(chainnan).  Welcoming  address  and  general  information. 

C.K.  Birdsall,  (2)T.L.  Crystal,  P.C.  Gray,  S.  Kuhn,  and  (l)Wm.S.  Lawson 

Introduction  to  PDWl  and  Q-machine  simulation;  presentation  of  “The  Movie” . 

C. K.  Birdsall,  (l)T.L.  Crystal,  P.C.  Gray,  P.  Knimm,  (2)S.  Kuhn,  Wm.S.  Lawson, 
M.  Oertl,  and  N.  Schupfer.  Trapped- electron  effects  on  negative-bias  d.c.  states  of  a 
collisionless  single -emitter  plasma  device:  theory,  simulation,  and  experiment. 

10:30  Coffee  and  posters  of  preceding  talks. 

11:15  K.S.  Theilhaber.  Dynamics  of  the  crossed-field  sheath:  a  movie  from  the  ES2  code. 

S.  Parker.  Electrostatic  potentials  due  to  nonuniform  magnetic  fields. 

D.  Hewett  (LLNL). 

Electromagnetic  boundary  conditions  for  implicit  codes. 

12:15  Lunch  in  Kerr  Campus  dining  hall  Building  10  (the  great  hall).  Posters. 

K.S.  Theilhaber  (chairman). 

1:30  (l)T.L.  Crystal,  M.  Horhager,  (2)S.  Kuhn,  and  (3)Wm.S.  Lawson. 

Theory  and  simulation  of  Pierce  diode  dynamics: 

Linear  and  nonlinear  features,  plus  nontrivial  external  circuit  effects. 

(l)C.K.  Birdsall  and  P.C.  Gray. 

Small- amplitude  impedance  Z(uj)  of  a  floating  single- emitter  plasma  device. 

2:30  Coffee  and  posters. 

3:30  C.K.  Birdsall,  B.I.  Cohen,  and  (1)R.  Procassini. 

Efficient  incorporation  of  Coulomb  collisions  in  bounded-plasma  simulations. 

L.A.  Schwager.  Transport  in  the  plasma  sheath  region  near  a  wall. 

S.  Kuhn.  Integral- equation  approaches  to  bounded- plasma  physics. 

5:00  Coffee  and  posters. 

6:30  Dinner  in  executive  dining  room  Building  10,  Kerr  Campus. 


Authors  listed  alphabetically 
(1),(2)  means  talks  split 


Schedule  of  BPPW  Presentations 
Tuesday  23  September  1986 


8:30  Coffee,  tea,  orange  juice  and  rolls. 

T.L.  Crystal(chairman). 

9:00  (1)  C.K.  Birdsall  and  K.S.  Th'eilhaber. 

Source-sheath  ion  acceleration 

P.C.  Gray,  (2)Wm.S.  Lawson  and  (l)R.D.  Pierce. 

I  —  V  curves  of  the  Child-Langmuir  electron  diode  from  particle  simulations. 

10:30  Coffee  and  posters  of  preceding  talks. 

(1)A.  Wendt,  M.A.  Lieberman  and  H.  Meuth. 

Theory  of  nonuniform  plasma  in  a  planar  magnetron  discharge 

(1)  J.  Helmer  and  K.  Doniger. 

Ion  diagnostics  in  a  planar  magnetron  discharge. 

12:15  Lunch. 

S.  Kuhn(chainnan). 

1:30  T.  Intrator  (Wisconsin). 

Some  comments  from  experiments. 

(l)Chris  Goedde,  M.A.  Lieberman,  A.J.  Lichtenberg. 

Stochastic  electron  heating  in  the  sheaths  of  an  RF  discharge. 

O.  Buneman  (Stanford). 

Matching  plasma  radiation  to  free-space  in  spherical  geometry  using  spherical 
harmonics  (outward  radiation  only). 

2:30  Coffee  and  posters. 

3:30  O.  Buneman  (Stanford). 

How  to  use  a  plasma  boundary  condition  stated  as  Z(u>)  in  plasma  simulation. 

Panel  discussion.  M.A.  Lieberman(chairman). 

Concluding  remarks  from  the  steering  committee. . .  short!  . . .  and  open  discussion. 
Where  are  we,  where  do  we  want  to  go,  and  what  should  we  do  next ? 

=> Should  we  have  another  BPPW  in  1988? 


SPECIAL  SECTION  V: 


Plasma  Theory  &  Simulation  Research  Group 
Prof.  C.K.  Birdsall 

EECS  Department  &  Electronics  Research  Laboratory 
University  of  California 
Berkeley,  CA  94720 


Publications  on 

BOUNDED-PLASMA  PHYSICS  AND  ENGINEERING 

Our  group  began  a  transition  into  bounded  plasma  problems  in  1981.  Our  new  objective  is  to 
understand  the  dynamics  of  plasma  sheaths  at  walls  in  terms  of  equilibrium,  stability  and  transport.  V*'e 
also  study  other  large  potentials  in  plasmas,  such  as  double  layers  and  thermal  barriers.  The  list  here 
reflects  the  transition,  a  learning  experience,  starting  with  reports  and  talks,  maturing  into  journal  arti¬ 
cles.  Major  suppon  comes  from  DOE  and  ONR,  plus  help  from  Varian  Associates-MICRO  and 
Hughes. 


1982 


ERL  Reports 

C.K.  Birdsall,  "Plasma-Sheath-Wall  Time-Dependent  Behavior  An  Informal  Survey,”  Memo.  No 
UCB/ERL  M82/51,  22  June  1982,  ERL,  UC  Berkeley.  (58  pages).  Talk  at  seminar  lead  by  Dr.  H. 
Ikegami  at  Biwa-ko,  Japan,  November  20,21,  1981. 


Conference  Proceedings 

C.K.  Birdsall,  W.S.  Lawson,  A.B.  Langdon,  “Plasma-surface  boundary  conditions;  plasma  diodes,” 
Proc.  of  Tenth  Conf.  on  Numerical  Simulation  of  Plasmas,  San  Diego,  CA,  Jan.  4-6,  1982. 

CJC.  Birdsall,  “Sheath  Formation  and  Fluctuations  with  Dynamic  Electrons  and  Ions,"  presented  at 
International  Conference  on  Plasma  Physics,  June  9-15,  1982,  Goteborg,  Sweden;  also  appears  in 
Conference  proceedings. 

CJC.  Birdsall,  "Electron  Diode  Dynamics;  Limiting  Currents;  Plasma  Diodes,”  presented  at  Symposium 
on  Plasma  Double  Layers,  Riso  National  Lab.,  Roskilde,  Denmark;  also  appears  in  Symposium  proceed¬ 
ings.  June  16-18.  1982. 

CJC.  Birdsall,  "Plasma  Sheath  Formation,  Ion  Acceleration,  and  Fluctuations  in  ‘Steady’  State,” 
presented  at  IEEE  International  Conference  on  Plasma  Science,  Ottawa,  Canada,  May  17-19,  1982. 
Paper  1R7 

Poster  PapeiVTalks 


Division  of  Plasma  Physics  Annual  Meeting,  Nov.  1982,  New  Orleans: 

CJC.  Birdsall,  “Temporal  development  of  a  plasma  sheath;  oscillatory  vs.  steady  state,"  Bull.  A.P.S. 
27,  page  1058,  October  1982. 

cf/updated  July  1987;  filename:  bounded. pubs 


1983 


Journal  Articles 

K.Y.  Kim,  “Weak  Monotonic  Double  Layers,"  Physics  Letters,  97  A,  8  August  1983. 

ERL  Reports 

Stephane  Rousset,  “Time-Dependent  Child-Langmuir  Diode  Simulation,"  Memo.  UCB/ERL  M83/39, 
11  July  1983. 

K.Y.  Kim,  “Theory  of  Asymmetric  Double  Layers,"  Memo.  UCB/ERL  M83/45,  22  July  1983. 

S.  Kuhn,  “Linear  Longitudinal  Oscillations  in  Collisionless  Plasma  Diodes  With  Thin  Sheaths.  Part  I. 
Method".  Memo.  UCB/ERL  M83/52,  31  August  1983. 

S.  Kuhn,  “Linear  Longitudinal  Oscillations  in  Collisionless  Plasma  Diodes  With  Thin  Sheaths.  Part  H 
Application  to  an  Extended  Pierce-Type  Problem,”  Memo.  UCB.ERL  M83/61,  3  October  1983. 

Poster  Papers/Talks 

Tenth  Conference  on  Numerical  Simulation  of  Plasmas  January  4-6,  1983,  San  Diego,  CA: 

(1)  C.K.  Birdsall,  W.S.  Lawson,  A.B.  Langdon  (LLNL),  “Plasma-Surface  Boundary  Conditions;Plasma 
Diodes.” 

(2)  W.S.  Lawson,  “A  Technique  for  Non-Periodic  Vlasov  Grid  Simulation.” 

APS/Division  of  Plasma  in  Physics,  November  1983,  Los  Angeles,  CA: 

(1)  S.  Kuhn,  “The  General  Perturbation  Problem  for  the  One-Dimensional  Collisionless  Plasma 
Diode." 

(2)  C.K.  Birdsall,  T.L.  Crystal,  S.Kuhn,  A.B.  Langdon,  W.S.  Lawson,  NP.  Otani,  “Plasma  Diode 
Simulation:  The  PDW1  Code.” 

(3)  Wm.  S.  Lawson,  “Vlasov  Simulation  on  a  Finite  Grid  With  A  Discontinuous  Distribution  Func¬ 
tion.” 

(4)  N.Otani,  “Scaling  of  Magnetized  Double  Layers.” 

(5)  K.Y.  Kim,  “Weak  Asymmetric  Double  Layers  (ADL).” 

(6)  TL.  Crystal,  S.  Kuhn,  Wm.  Lawson,  "Pierce  Diode  Instabilities  and  Extensions:  Simulations  With 
PDW1.” 


Journal  Articles 


S.Kuhn,  “Linear  Longitudinal  Oscillations  in  Collisionless  Diodes  with  Thm  Sheaths;  Pan  I.  Method," 
Physics  of  Fluids  27(7),  1821,  July. 

S.Kuhn.  “Linear  Longitudinal  Oscillations  in  Collisionless  Plasma  Diodes  with  Thin  Sheaths;  Pan  II. 
Application  to  an  Extended  Pierce-type  Problem,”  Physics  of  Fluids27(l),  1834,  July. 


ERL  Reports 

W.S.  Lawson,  “A  Method  for  Simulation  of  Vlasov  Systems  with  Discontinuous  Distnbution  Functions 
and  GASBAG  User's  Manual,”  Memo  No.  UCB/ERL  M84/24,  March,  UC  Berkeley. 

W.S  Lawson,  “PWD1  User's  Manual,"  Memo  No.  UCB/ERL  M84/37,  Apnl,  UC  Berkeley. 

K.Y.  Kim,  “ Vlasov-Poisson  and  Modified  k-dV  Theory  and  Simulation  of  Weak  and  Strong  Double 
Layers,"  Memo  No.  UCB.ERL  MS4/47,  June,  UC  Berkeley. 

T.L.  Crystal  and  S.  Kuhn,  “Particle  simulation  of  the  low -alpha  Pierce  diode,"  Memo  No.  UCB/ERL 
MS4/74,  20  September  1984,  UC  Berkeley. 


Poster  Papers/Talks 


International  Conference  on  Plasma  Physics,  June  27-July  3,  1984,  Lausanne,  Switzerland: 

(1)  S.  Kuhn,  C.K.  Birdsall,  TL.  Crystal,  Pi-.  Gray,  William  S.  Lawson,  “One  dimensional  particle 
simulations  of  the  collisionless  single-ended  Q  machine,"  (P8-11). 

(2)  TL..  Crystal  and  S.  Kuhn,  "Particle  Simulations  of  the  low-alpha  Pierce  diode,"  (P33-11). 


Second  Symposium  on  Plasma  Double  Layers  and  Related  Topics,  July  5-6,  1984,  Innsbruck,  Austria: 

(1)  P.  Gray,  S.  Kuhn,  TL.  Crystal  and  C.K.  Birdsall,  “Particle  simulations  of  the  unstable  states  of  a 
collisionless  single-end  plasma  device." 

(2)  NT.  Otani,  "Evolution  of  ion  holes  and  large  potential  structures  in  simulations  of  bounded 
cunem-driven  magnetized  plasma  systems.” 

(3)  K.Y.  Kim  and  TL.  Crystal,  "Theory  and  Simulation  of  Ion  Acoustic  Double  Layers." 


APS  Division  of  Plasma  Physics  Meeting,  November  1984,  Boston,  MA: 

(1)  C.K.  Birdsail,  William  S.  Lawson,  T.L.  Crystal,  S.  Kuhn,  A.B.  Langdon,  “Plasma  diode  simulation: 
the  PDW1  code.” 

(2)  L.A.  Schwager.  “Simulation  of  charged  particle  reflection  and  emission  effects  in  the  plasma- 
sheath  region.” 

IEEE  Conference  on  Plasma  Science,  May  1984,  St.  Louis,  MO: 

(1)  S.  Kuhn,  P.Gray,  TJ-.  Crystal,  C.K.  Birdsail,  “Axial  Equilibria  of  Collisionless  Single-Ended  Q 
Machines:  Panicle  Simulations  versus  Theory  and  Experiment.”  (10  minute  movie) 

(2)  TJL.  Crystal,  S.  Kuhn  and  C.K,  Birdsail,  “The  Classical  Pierce  Diode:  Using  Particle  Simulations 
on  Linear  and  Nonlinear  Behavior  and  Final  States.” 
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1985 


Book 

CJC.  Birdsall,  A.B.  Langdon,  Plasma  Physics  via  Computer  Simulation,  McGraw-Hill,  New  York.  See 
especially  Chapters  14,  15,  16  on  bounded  plasmas. 

Journal  Articles 

W.S.  Lawson,  “A  Method  for  Simulation  of  Vlasov  Systems  with  Discontinuous  Distribution  Func¬ 
tions,”  J.  of  Computational  Physics,  61  (1),  p.  51,  October. 

S.  Ishiguro  and  T.  Kamimura,  ‘‘Double  Layer  Formation  Caused  by  Contact  Between  Different  Tem¬ 
perature  Plasmas,"  Physics  of  Fluids  28  (7),  pp.  2100-2105,  July. 

T. L.  Crystal  and  S.  Kuhn,  ‘‘Particle  Simulations  of  the  Low-alpha  Pierce  Diode,”  Physics  of  Fluids  28 
(7),  pp.  2116-2124,  July. 


Reports 


Bruce  I.  Cohen  and  Mark  E.  Stewart;  Charles  K.  Birdsall,  "Direct  Implicit  Particle  Simulation  of  Tan¬ 
dem  Mirrors, "  LLN’L  Mirror  Theory  Monthly  February  1985. 

C.  K.  Birdsall,  T.  L.  Crystal,  P.  C.  Gray  and  S.  Kuhn,  "Bounded  Plasma  Dynamics  from  Particle  Simu¬ 
lations;  Movie  Script,"  Memo  No.  UCB/ERL  M85/54,  2  July  1985,  UC  Berkeley. 


Poster  Papers/Talks 

Sherwood  Theory  Conference,  April  15-17,  1985,  University  of  Wisconsin,  Madison,  Wisconsin: 

L.A.  Schwager  and  C.K.  Birdsall,  “A  Model  of  the  Plasma-Sheath  Region  Including  Secondary  Elec¬ 
tron  Emission  and  Ion  Reflection.” 


Eleventh  International  Conference  on  Numerical  Simulation  of  Plasma,  Montreal,  Quebec,  Canada,  25- 
27  June  1985: 

(1)  CJC.  Birdsall,  W.S.  Lawson,  T.  Crystal,  S.  Kuhn,  N.  Otani,  I.  Roth,  and  A.B.  Langdon,  “Problems 
of  Bounded  Particle  Simulations  and  First  Generation  Solutions.” 

(2)  B.I.  Cohen,  M.E.  Stewart,  and  C.K.  Birdsall,  “Direct  Implicit  Particle  Simulation  of  Tandem  Mir¬ 
ror.” 


17 ih  International  Conference  on  Phenomena  in  Ionized  Gases,  Budapest,  July  1985: 

CJC.  Birdsall,  Ti_.  Crystal,  P,  Gray,  and  S.  Kuhn,  “Instability  Of  A  Collisionless  Single-Ended  Plasma 
Device  Operated  In  The  Positive-Bias  Electron-Rich  Regime;  From  Simulations." 


INTOR  Specialists'  Meeting  On  Impurity  Control,  IAEA  Headquarters,  Vienna,  Austria; 
C.K.  Birdsall,  “Boundary  Conditions;  Matching  Plasma-Sheath-Wall.”  (invited  paper) 


APS/Division  of  Plasma  Physics  Annual  Meeting,  Nov  4-8,  1985,  San  Diego,  CA: 

(1)  B.I.  Cohen,  M£.  Stewart  RP.  Freis  (LLNL),  L.A.  Strugala,  and  C.K.  Birdsall,  “Direct  Implicit 
Particle  Simulation  of  Tandem  Mirrors.” 

(2)  Wm.  Lawson,  P.C.  Gray,  EA.  Adler,  and  C.K.  Birdsall,  “A  Simulation  Study  of  Effects  of  Emis¬ 
sion  Noise  in  a  Thermionic  Diode." 

(3)  P.C.  Gray,  T.L.  Crystal,  S.  Kuhn,  and  C.K.  Birdsall,  “  Equilibrium  Corrections  for  Electron  Hole 
Filling:  Simulation  vs.  Theory." 

(4)  L.A.  Schwager  and  C.K.  Birdsall,  "A  Model  of  the  Plasma-Sheath  Region  Including  Secondary 
Electron  Emission  and  Ion  Reflection.” 

(5)  A.  Wendt,  P.C.  Gray,  H.  Meuth,  M.A.  Lieberman  and  C.K.  Birdsall,  ‘..‘Experimental  and  Simulation 
Study  of  Planar  Magnetron  Discharges." 

(6)  C.K.  Birdsall,  P.C.  Gray  and  T.L.  Crystal,  "A  Simulation  Study  of  the  Dependence  of  the  Plasma- 
Sheath  Impedance  on  Frequency." 


1986 


Journal  Articles 

Charles  K.  Birdsall,  Niels  F.  Otani  and  Bruce  1.  Cohen,  "Simulation  of  Plasma  Dynamics  Using  Many 
Particles,"  IEEE  Electrotechnology  Review  2,  pp.  43,44,  1986. 

S.  Kuhn  and  M.  HSrhager,  “External-circuit  Effects  on  Pierce-Diode  Stability  Behavior,”  J.  Appl.  Phys. 
60  (6),  pp.  1952-1959,  September  15,  1986. 


Reports 


William  S.  Lawson,  "Linear  Magnetized  Plasma  Response  to  an  Oblique  Electrostatic  Wave,"  Univer¬ 
sity  of  California,  Berkeley,  Memorandum  No.  UCB/ERL  M86/94,  December  15,  1986. 


Poster  Papers 

APS  Division  of  Plasma  Physics  Twenty-Eighth  Annual  Meeting,  November  3-7,  1986,  Baltimore,  MD: 

L.A.  Schwager  and  C.K.  Birdsall,  “Transport  across  the  Plasma-Wall  Sheath  with  Secondary  Electron 
Emission  and  Ion  Reflection." 

T.L.  Crystal  and  S.  Kuhn,  “  Pierce  Diode  Instability  Simulations  with  External  Inductance  and  Resis¬ 
tance." 

S-E.  Parker,  C.K.  Birdsall,  and  K.  Theilhaber,  “Electrostatic  Effects  on  Confinement  Outside  the 
Separatrix  of  Field  Reversed  Configurations.” 

RJ.  Procassini,  C.K.  Birdsall,  and  B.I.  Cohen  (LLNL),  “Direct  Implicit  Particle  Simulation  of  Tandem 
Mirrors.” 

K.  Theilhaber  and  C.K.  Birdsall,  “Structure  of  the  Qossed-Field  Sheath.”  (with  short  movie) 


C.  K.  Birdsall  and  the  Plasma  Theory  and  Simulation  Group,  "Planar  Plasma  Discharges:  Simulations 
and  Experiments,"  at  Sandia  National  Labs,  October  7,  1986. 

C.  K.  Birdsall  and  Michael  Lieberman,  "Planar  Plasma  Discharges:  Simulations  and  Experiments”  at 
Varian  Associates  Technical  Seminar,  October  20,  1986. 


Workshop 

Bounded-Plasma  Physics  Workshop,  University  of  California,  Berkeley,  September  22-23,  1986.  See 
schedule  following. 


Bounded-Plasma  Physics  Workshop 
University  of  California,  Berkeley 
September  22-23,  1986 


Schedule  of  BPPW  Presentations 
Monday  22  September  1986 

8:30  Registration  and  parking  arrangements. 

Coffee,  tea,  orange  juice  and  rolls. 

9:00  C.K.  Birdsall(chairman).  Welcoming  address  and  general  information. 

C.K.  Birdsall,  (2)T.L.  Crystal,  P.C.  Gray,  S.  Kuhn,  and  (l)Wm.S.  Lawson 

Introduction  to  PDWl  and  Q-machine  simulation;  presentation  of  “The  Movie". 

C.K.  Birdsall,  (l)T.L.  Crystal,  P.C.  Gray,  P.  Krnmm,  (2)S.  Kuhn,  Wm.S.  Lawson, 
M.  Oertl,  and  N*.  Schupfer.  Trapped- electron  effects  on  negative-bias  d.c.  states  of  a 
colhsionless  single -emitter  plasma  device:  theory,  simulation,  and  experiment. 

10:30  Coffee  and  posters  of  preceding  talks. 


11:15  K.S.  Theilhaber.  Dynamics  of  the  crossed-field  sheath:  a  movie  from  the  ES2  code. 
S.  Parker.  Electrostatic  potentials  due  to  nonuniform  magnetic  fields. 

D.  Hewett  (LLNL). 

Electromagnetic  boundary  conditions  for  implicit  codes. 

12:15  Lunch  in  Kerr  Campus  dining  hall  Building  10  (the  great  hall).  Posters. 

K.S.  Theilhaber(chairman). 

1:30  (l)T.L.  Crystal,  M.  Horhager,  (2)S.  Kuhn,  and  (3)Wm.S.  Lawson. 

Theory  and  simulation  of  Pierce  diode  dynamics: 

Linear  and  nonlinear  features,  plus  nontrivial  external  circuit  effects. 

(l)C.K.  Birdsall  and  P.C.  Gray. 

Small- amplitude  impedance  Z(-j)  of  a  floating  single- emitter  plasma  device. 
2:30  Coffee  and  posters. 

3:30  C.K.  Birdsall,  B.I.  Cohen,  and  (1)R.  Procassini. 

Efficient  incorporation  of  Coulomb  collisions  in  bounded- plasma  simulations. 

L.A.  Schwager.  Transport  »n  the  plasma  sheath  region  near  a  wall. 

S.  Kuhn.  Integral- equation  approaches  to  bounded- plasma  physics. 

5:00  Coffee  and  posters. 

6:30  Dinner  in  executive  dining  room  Building  10,  Kerr  Campus. 


Authors  listed  alohabetically 
(1),(2)  means  talks  split 
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Schedule  of  BPPW  Presentations 
Tuesday  23  September  1986 


8:30  Coffee,  tea,  orange  juice  and  rolls. 

T.L.  Crystal  (chairman). 

9:00  (1)  C.K  Birdsall  and  K.S.  Theilhaber. 

Source-sheath  ion  acceleration 

P.C.  Gray,  (2)\Vm.S.  Law-son  and  (l)R.D.  Pierce. 

I  —  V  curves  of  the  Child- Langmuir  electron  diode  from  particle  simulations. 


10:30  Coffee  and  posters  of  preceding  talks. 


(1)A.  Wendt,  M.A.  Liebennan  and  H.  Meuth. 

Theory  of  nonuniform  plasma  in  a  planar  magnetron  discharge 

(1)  J.  Helmer  and  K.  Doniger. 

Ion  diagnostics  in  a  planar  magnetron  discharge. 


12:15  Lunch. 

S.  Kuhn  (chairman). 

1:30  T.  Intrator  (Wisconsin). 

Some  comments  from  experiments. 

(l)Chris  Goedde,  M.A.  Liebennan,  A.J.  Lichtenberg. 

Stochastic  electron  heating  in  the  sheaths  of  an  RF  discharge. 

O.  Buneman  (Stanford). 

Matching  plasma  radiation  to  free-space  in  spherical  geometry  using  spherical 
harmonics  (outward  radiation  only). 

2.30  Coffee  and  posters. 

3:30  O.  Buneman  (Stamford). 

How  to  use  a  plasma  boundary  condition  stated  as  Z(lj)  in  plasma  simulation. 

Panel  discussion.  M.A.  Lieberm  an  (chairman). 

Concluding  remarks  from  the  steering  committee. . .  short!  . . .  and  open  discussion. 
Where  are  we,  where  do  we  want  to  go,  and  what  should  we  do  next? 

=>Should  we  have  another  BPPW  in  1988? 


1987 


Reports 

William  S.  Lawson,  "Computer  Simulation  of  Bounded  Plasma  Systems,"  University  of  California, 
Berkeley,  Memorandum  No.  UCB/ERL  M87/14,  March  5,  1987. 

K.  Theilhaber  and  C.  K.  Birdsall,  "Vortex  Dynamics  and  Transport  to  the  Wall  in  a  Crossed-Field 
Plasma  Sheath,"  University  of  California,  Berkeley,  Memorandum  No.  UCB/ERL  M87/18,  Apnl  10, 
1987. 

Kim  Theilhaber,  "ES2  User’s  Manual-Version  1,"  University  of  California,  Berkeley,  Memorandum 
No.  UCB/ERL  MS '.'23,  May  1  1,  1987. 

R.  J.  Procassini,  C.  K.  Birdsall,  E.  C.  Morse,  and  B.  I.  Cohen,  "A  Relativistic  Monte  Carlo  Binary  Col¬ 
lision  Model  for  Use  in  Plasma  Particle  Simulation  Codes,"  University  of  California,  Berkeley, 
Memorandum  No.  UCB/ERL  M87/24,  May  14, 1987. 

William  S.  Lawson,  "Artificial  Cooling  Due  to  Quiet  Injection  in  Bounded  Plasma  Particle  Simula¬ 
tions,"  University  of  California,  Berkeley,  Memorandum  No.  UCB/ERL  M87/34,  May  21,  1987. 

William  S.  Lawson,  "The  Pierce  Diode  with  an  External  Circuit.  I.  Simulation  in  the  Linear  Regime," 
University  of  California,  Berkeley,  Memorandum  No.  UCB/ERL  MS7/51,  July  1987. 

William  S.  Lawson,  "The  Pierce  Diode  with  an  External  Circuit.  II.  Non-uniform  Equilibria,"  Univer¬ 
sity  of  California,  Berkeley,  Memorandum  No.  UCB/ERL  M87/52,  July  1987. 

William  S.  Lawson,  "The  Pierce  Diode  with  an  External  Circuit.  III.  Chaotic  Behavior,”  to  be  ERL 
Memo.,  1987. 

Lou  Ann  Schwager,  "Effects  of  Secondary  Electron  Emission  on  the  Collector  and  Source  Sheaths 
Bounding  a  Finite  Ion  Temperature  Plasma"  to  be  ERL  Memo.,  1987. 

Lou  Ann  Schwager,  "Effects  of  Ion  Reflection  on  the  Collector  and  Source  Sheaths  Bounding  a  Finite 
Ion  Temperature  Plasma"  to  be  ERL  Memo.,  1987. 


Poster  Papers/Talks 

Poster  Paper  ai  US. /Japan  Seminar:  Effects  of  Electric  Fields  on  Magnetic  Confinement.  January  22- 
24,  1987,  University  of  California  Los  Angeles: 

K.  S.  Theilhaber  and  C.  K.  Birdsall,  "Large  Electric  Fields  in  a  Magnetized  Plasma  Sheath:  Long-lived 
Vortices." 

Poster  Papers  at  Sherwood  Controlled  Fusion  Theory  Conference.  April  6-8,  1987,  San  Diego,  Califor¬ 
nia: 

K.  Theilhaber,  "Transport  Induced  by  a  Crossed-Field  Sheath." 

Richard  J.  Procassini  and  Charles  K.  Birdsall,  "Performance  and  Optimization  of  Direct  Implicit  Time 
Integration  Schemes  for  Use  in  Electrostatic  Particle  Simulation  Codes." 
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Lou  Ann  Schwager,  "Collector  Sheath  and  Source  Sheath  in  a  Collisionless  Finite  Ion  Temperature 
Plasma  with  Secondary  Electron  Emission  and  Ion  Reflection  at  the  Bounding  Surface." 

Poster  Papers  at  IEEE  International  Conference  on  Plasma  Science,  June  1-3,  1987,  Arlington,  Va.: 

K.  Theilhaber  and  C.  K.  Birdsall,  "Large  Electric  Fields  in  a  Magnetized  Plasma  Sheath;  Long-lived 
Vortices." 

L.  A.  Schwager  and  C.  K.  Birdsall,  "Potential  Drop  and  Transport  in  a  Bounded  Plasma  Including 
Secondary  Electron  Emission  and  Ion  Reflection  at  the  Collector." 

ln\ited  talk  at  APS  Division  of  Plasma  Physics  Twenty-Ninth  Annual  Meeting,  November  2-6,  1987, 
San  Diego,  California: 

K.  Theilhaber,  "Vortex  Formation  and  Transport  to  the  Wall  in  a  Cxossed-Field  Sheath." 

Poster  Papers  at  APS  Division  of  Plasma  Physics  Twenty-Ninth  Annual  Meeting,  November  2-6,  1987, 
San  Diego,  California: 

C.  K.  Birdsall,  K.  S.  Theilhaber,  and  S.  Kuhn,  "Ion  Acceleration  in  a  Plasma  Source  Sheath." 

A.  Friedman  and  S.  L.  Ray,  LLNL,  and  S.  Parker  and  C.  K.  Birdsall.  U.  C.  Berkeley,  "Prospects  for 
Multi-Scale  Particle-in-Cell  Simulation  of  Plasmas." 

Wm.  S.  Lawson,  "Investigations  of  the  Pierce  Diode  Strange  Attractor." 

S.  E.  Parker  and  C.K.  Birdsall,  V.  C.  Berkeley,  and  A.  Friedman  and  S.  L.  Ray,  LLNL,  "Direct  Implicit 
Particle  Simulation  of  a  Bounded  Plasma  System." 

R.  J.  Procassini,  J.  C.  Cummings  and  C.  K.  Birdsall,  University  of  California,  Berkeley,  and  B.  I. 
Cohen,  Lawrence  Livermore  National  Laboratory,  "Direct  Implicit  Particle  Simulation  of  Simple  Mir¬ 
rors." 

L.  A.  Schwager.  "The  Effect  of  Thermionic  and  Secondary  Electron  Emission  at  the  Collector  on  Poten¬ 
tial  Drop  and  Transport  Through  the  Plasma-Sheath  Region." 
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